ChinaXiv 合 作 期 刊 


202302.00194v1 


chinaXiv 


时 变 多 车 型 下 的 生 鲜 农产品 配送 路 径 优化 模型 


EA, RAA”, R 


Rak MaE R 


(1. 广西 大 学 计算 机 与 电子 信息 学 院 ,广西 南宁 $30004; 2. 国家 农业 信息 化 工程 技术 研究 中 心 ,北京 100097) 


摘 要 : 针对 农业 供应 链 的 运输 环节 中 生 鲜 农产品 配送 模型 存在 的 速度 恒定 、 碳 排放 计算 方法 单一 的 问题 ， 
本 研究 结合 路 网 时 变 特征 和 新 的 多 车 型 碳 排 放 计算 方法 ， 提 出 了 考虑 配送 距离 、 多 车 型 碳 排 放量 、 货 物 损 
耗 和 和 车 辆 固定 成 本 等 4 个 优化 目标 的 生 鲜 农产品 配送 路 径 优化 模型 ;并 根据 模型 特点 提出 了 一 种 改进 的 双 策 


略 种 群 协同 蚁 群 算法 (Double-Strategies Co-Evolutionary Ant Colony System, DC-ACS )。 利 用 改进 蚁 群 算法 


对 Solomon 数据 集 的 C105 算 例 进行 了 求解 ， 在 4 个 优化 目标 上 分 别 取得 最 优 解 为 937.94 km、4961.48 元 、 


4081.78 元 和 7500.87 元 ,证 明了 本 研究 提出 的 模型 的 有 效 性 。 在 模型 有 效 的 基础 上 ， 通 过 试验 结果 证 明 ， 改 


进 蚁 群 算法 比 基 本 蚁 群 算法 在 4 个 优化 目标 上 的 配送 总 成 本 平均 降低 幅度 超过 14%， 证 明 改 进 蚁 群 算法 更 具 
有 优越 性 。 使 用 改进 疏 群 算法 对 集中 、 随 机 和 混合 3 种 不 同 分 布 的 大 规模 算 例 进行 求解 ，3 种 分 布 上 分 别 求 
得 最 优 总 成 本 为 19,939.53 、24,095.00 和 24,397.58 元 。 综 上 所 述 ， 所 提 模 型 和 算法 可 以 为 冷 链 物流 企业 的 城 
市 配送 路 径 决 策 提 供 良 好 的 参考 依据 ， 对 完善 智慧 农业 供应 链 的 配送 路 径 优化 模型 和 优化 方法 提供 了 新 的 


思路 ， 为 企业 进一步 扩大 规模 提供 了 参考 。 
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abstract) 


1 引言 


随 着 生活 水 平 的 提高 ， 人 们 对 鲜 肉 、 水 产品 
和 有 果蔬 等 生 鲜 农产品 的 需求 不 断 增加 ， 对 生 鲜 易 
腐 品 的 冷 链 配送 服务 的 要 求 日 益 增 加 。 但 城市 冷 
链 配送 易 受 天 气 和 交通 状况 的 影响 ， 进 而 导致 配 
送 成 本 和 运输 时 间 增 加 。 准 确 量化 道路 状况 对 速 
度 的 影响 是 城市 配送 模型 建立 中 需要 解决 的 问 
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题 。 而 客户 订单 数量 的 增多 和 分 布 的 复杂 化 " 
使 得 冷 链 企 业 的 运营 成 本 进一步 增加 ， 过 去 单一 
车 型 的 配送 模式 无 法 满足 企业 的 要 求 。 多 车 型 配 
送 模式 的 出 现 解决 了 这 一 问题 ”“， 计 算 配 送 过 程 
中 不 同类 型 车 辆 的 碳 排放 并 最 大 化 降低 碳 排 放 和 
企业 成 本 成 为 研究 的 热点 。 对 于 冷 链 物流 企业 来 
说 ,无 法 准确 量化 复杂 道路 状况 和 多 车 型 矶 排放 
对 配送 成 本 的 影响 ,制约 了 企业 进一步 发 展 ， 同 
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时 也 影响 高 效 、 绿 色 的 智 意 农业 供应 链 的 
建设 5。 

针对 模型 中 速度 固定 和 车 辆 碳 排 放 计 算 方法 
单一 问题 ， 国 内 外 许多 研究 者 进行 了 深入 人 研究 。 
陈 久 梅 等 “ 考虑 了 不 同 产 品 对 冷冻 温度 和 湿度 
的 不 同 需求 ， 建 立 了 生 鲜 农产品 多 陋室 车 辆 路 径 
优化 模型 ( Multi-Compartment Vehicle Routing 
Problem, MCVRP), ， 接 着 将 车 辆 路 径 规 划 问 题 
(Vehicle Routing Problem, VRP) 的 国际 标准 算 
例 改 变 后 得 到 28 个 MCVRP 标 准 算 例 ， 最 后 使 用 
融合 了 邻 域 搜索 算法 的 粒子 群 算 法 对 模型 求解 ， 
改进 后 的 算法 求 得 的 总 路 径 降 低 了 2%。 张 倩 
等 "建立 了 考虑 配送 成 本 、 产 品 新 鲜 度 、 碳 排 
放 等 因素 的 生 鲜 电 商 优化 模型 ， 其 中 计算 产品 新 
鲜 度 时 考虑 了 距离 和 时 间 因 素 形 成 新 鲜 度 公式 ， 
在 计算 碳 排放 时 考虑 了 运输 距离 和 负载 ， 构 建 了 
包括 总 成 本 、 产 品 新 鲜 度 和 碳 排放 的 多 目标 优化 
模型 ， 最 后 使 用 果 蝇 算法 对 该 模型 求解 。 上 述 研 
究 虽 然 构建 了 考虑 碳 排放 的 模型 ， 但 忽视 了 车 辆 
速度 变化 的 现实 情况 ， 同 时 在 碳 排放 计算 方法 上 
没有 创新 。 肖 建华 等 ”提出 了 基于 城市 道路 限 
行 的 多 车 型 车 辆 路 径 优化 模型 ， 并 使 用 变 邻 域 搜 
索 算 法 求解 该 模型 ， 证 明 该 模型 可 以 降低 物流 成 
本 、 提 升 配送 效率 和 改善 城市 环境 。 赵 志学 等 ” 
提出 了 考虑 拥堵 区 域 的 多 车 型 绿色 车 辆 路 径 问题 
优化 模型 ， 求 解 结果 表明 不 同 的 道路 状况 对 物流 
配送 成 本 有 和 较 大 影响 。 李 军 涛 等 " 提出 了 考虑 
拥堵 指数 的 多 车 型 冷 链 路 径 优 化 模型 ， 并 使 用 自 
适应 遗传 模拟 退货 算法 求解 ， 结 果 表明 该 模型 和 
算法 可 以 降低 56% 的 碳 排 放量 。Wang 等 "在 
2016 年 构建 了 总 成 本 最 低 和 最 大 新 鲜 度 的 多 目标 
优化 模型 ， 并 使 用 两 阶段 启发 式 算法 求解 ， 其 中 
先 使 用 K-means 聚 类 方法 生成 初始 解 ， 接 着 使 用 
可 变 邻 域 搜索 和 遗传 算法 对 其 进行 进一步 寻 优 。 
Ma 等 ”在 时 间 依 赖 车 辆 路 径 问 题 的 基础 上 
(Time Dependent Vehicle Routing Problem, TD- 
VRP)， 建 立 了 车 辆 灵活 出 发 条 件 下 的 生 鲜 农 产 
品 配送 模型 。Liu 等 ' 引 在 带 时 间 窗 的 时 间 依 赖 车 


辆 路 径 问 题 (Time Dependent Vehicle Routing 
Problem With Time Window, TDVRPTW) 的 基 
hE, Beth FAR TE AS. AJRAM 
耗 的 配送 路 径 优化 模型 ， 并 使 用 改进 的 蚁 群 算法 
对 该 模型 求解 。Hsiao 等 "利用 改进 遗传 算法 人 研 
究 了 视频 类 型 和 温度 变化 的 路 径 优化 问题 。 
Byung 等 “” 比较 了 冷藏 车 和 一 般 车 型 的 配送 成 
本 ,确认 了 冷藏 车 在 生 鲜 农产品 配送 中 的 有 效 
性 。Alinaghian 和 Nadia 等 ' 采用 大 规模 邻 域 算 
法 和 变 邻 域 算 法 解决 了 多 隔 室 和 多 中 心 的 配送 路 
径 问 题 。Hiassat 等 ” 考虑 易 腐 品 特性 ， 建 立 了 
配送 中 心 -库存 优化 模型 ， 并 使 用 局 部 搜索 改进 
的 遗传 算法 求解 。 以 上 研究 同时 考虑 了 多 车 型 和 
道路 拥堵 情况 ,但 对 于 晴天 时 不 同时 间 段 的 道路 
状况 ,没有 建立 起 相应 的 车 速 特征 模型 。 最 后 ， 
部 分 研究 1!“ 的 算 例 客户 点 数量 均 在 30 个 以 下 ， 
日 算 例 客户 分 布 类 型 较为 单一 ， 没 有 考虑 大 规模 
客户 在 离散 、 集 中 与 混合 分 布 的 这 三 种 情况 下 的 
求解 结果 。 

综 上 所 述 ， 当 前 有 关 冷 链 配 送 路 径 优化 的 研 
究 中 ， 在 模型 构建 上 ， 主 要 探讨 了 道路 拥堵 条 件 
下 的 模型 构建 和 求解 ， 没 有 考虑 到 现实 环境 中 车 
辆 速度 受 时 间 影 响 而 变化 的 实际 情况 “站 ， 冷 链 
配送 车 辆 碳 排 放 计 算 方 法 也 较为 单一 ; 在 模型 求 
解 中 ， 对 大 规模 和 不 同 客户 分 布 类 型 算 例 的 求解 
较 少 221。 基于 上 述 背 景 ， 本 人 研究 在 传统 的 生 鲜 
农产品 配送 路 径 优化 模型 的 基础 上 ， 考 虑 了 路 网 
时 变 特征 并 采用 了 基于 载重 的 多 车 型 碳 排放 计算 
模型 ， 提 出 综合 考虑 碳 排放 、 和 车 辆 路 径 、 货 物 新 
鲜 度 和 车 辆 固定 成 本 的 多 目标 优化 的 车 辆 路 径 模 
型 ; 通过 设计 改进 蚁 群 算法 来 求解 大 规模 生 侠 农 
产品 配送 算 例 。 为 冷 链 物流 企业 的 城市 配送 路 径 
决策 提供 良好 的 参考 依据 。 


2 问题 描述 与 模型 构建 


2.1 问题 描述 和 假设 
本 研究 探索 车 辆 速度 随时 间 变 化 条 件 下 的 多 
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车 型 物流 配送 车 辆 路 径 问 题 ， 具 体 描 述 为 : 物流 


时 间 段 随 着 逐渐 接近 早晚 高 峰 ， 和 车辆 受 影响 程度 


配送 中 心 具 有 多 种 类 型 的 车 辆 ， 给 位 于 城市 不 同 
地 点 的 客户 配送 货物 ; 客户 点 的 货物 需求 量 和 位 
置 确定 ， 具 有 时 间 窗 要 求 ; 随 着 时 间 的 变化 ， 道 
路 状况 会 影响 车 辆 的 行驶 速度 。 优 化 目标 是 在 考 
虑 车 辆 油耗 、 碳 排放 、 和 车辆 固定 成 本 和 货 损 成 本 
的 情况 下 ， 找 出 满足 客户 特定 要 求 的 车 辆 调度 和 
路 径 安 排 方案 。 

为 方便 对 问题 分 析 和 求解 ， 做 出 以 下 假设 : 

(1) 配送 中 心 和 客户 点 的 位 置 坐标 已 知 ， 车 
辆 的 起 点 和 终点 均 为 配送 中 心 。 

(2) 配送 中 心 有 不 同类 型 的 配送 车 辆 ， 且 不 
同类 型 车 辆 的 数量 有 限 。 

(3) 车 辆 行驶 在 不 同 的 时 间 段 ， 速 度 不 同 。 

(4) 每 个 客户 配送 时 间 段 固定 ， 提 前 或 者 超 
过 服务 时 间 窗 均 要 承担 惩罚 费用 。 

(5) 每 个 客户 由 一 辆 车 为 其 提供 服务 ， 客 户 
的 产品 需求 量 确定 ， 均 小 于 配送 车 辆 的 最 大 承载 
能 力 。 

(6) 不 同类 型 车 辆 的 最 大 承载 能 力 、 固 定 成 
本 费用 、 油 耗 和 碳 排放 率 不 同 。 

(7) 配送 车 辆 在 等 待 和 服务 期 间 ， 没 有 油耗 
和 碳 排放 产生 。 


2.2 模型 构建 


2.2.1 车 辆 行驶 时 间 分 析 

配送 车 辆 的 速度 受到 较 多 因素 的 影响 ， 在 传 
统 的 车 辆 路 径 优 化 配送 模型 中 ， 道 路 状况 和 天 气 
保持 不 变 ， 进 而 配送 车 辆 速度 恒定 ; 在 实际 的 配 
送 过 程 中 ， 变 化 的 行驶 速度 更 符合 实际 。 因 此 ， 
本 研究 分 析 了 在 晴天 状况 时 ， 早 晚 高 峰 对 车 辆 行 
驶 速度 的 影响 。 

文献 [23] 建立 了 考虑 实时 路 况 和 接 驶 点 的 
冷 链 物流 配送 路 径 优 化 模型 ， 其 中 的 实时 路 况 采 
用 道路 实时 路 况 系数 表述 ; 文献 [21,24] 根据 
天 气 状 况 的 不 同 ， 建 立 了 车 速 特征 模型 。 本 研究 
参考 上 述 研究 成 果 和 相关 监测 数据 后， 建立 
了 晴天 时 与 不 同时 间 段 下 的 车 速 特征 模型 ， 如 公 
式 (1) 所 示 。 在 每 天 6:00 一 8:00 和 12:00 一 18:00 


不 断 增 大 ; 在 8:00 一 12:00 和 18:00 一 22:00 时 间 
段 ， 道 路 上 车 辆 数 减少 的 同时 ， 对 配送 车 辆 的 影 
响 也 逐渐 降低 ; 其 他 时 间 段 无 影响 。 

青天 时 不 同时 段 对 车 辆 速度 的 影响 率 , 表 


IRA: 
Can = 
0 0 < t < 6 或 22 < 1<24 
0.225t - 1.35 6<t<8 
-0.05¢ + 0.85 8<t<12 
t 
—-0.15 12<t<18 
30 
-—0.1125t¢ + 2.475 18<t<22 


(1) 
其 中 ,1 为 车 辆 行驶 时 所 在 的 时 间 点 。 
青天 时 计算 不 同时 段 下 的 车 辆 速 度 $ 如 
公式 (2): 


S,=8,-(-£,) (2) 
其 中 ，5 ,是 配送 车 辆 的 平均 行驶 速度 ，km/h。 
因此 ， 配 送 车 辆 在 从 任意 点 i 到 j 的 行驶 时 间 表 
示 为 : 
ti aii (3) 
Š, 


其 中 ， 力 为 配送 车 辆 从 ;点 到 /点 需要 花费 的 
时 间 ，h; 4 为 任意 客户 点 i 到 j 的 距离 ，km。 
2.2.2 ”考虑 载重 的 油耗 与 碳 排 放 计 算 

不 同 车 型 的 营运 车 辆 ， 在 计算 油耗 和 碳 排 放 
时 会 受到 多 种 因素 的 影响 ， 其 中 载重 是 最 重要 的 
影响 因素 。 首 先 经 过 测算 得 到 不 同 载重 货车 的 百 
公里 油耗 均值 ， 每 种 类 型 货车 取 不 同 数量 ， 具 体 
见 表 1。 

接着 根据 碳 平 衡 原 理 计算 各 类 货车 的 CO, 排 
WAT, 最 后 取 各 货车 的 结果 均值 作为 车 辆 的 

表 1 试验 车 型 信息 


Table 1 Information of the test vehicles 


车 型 额定 载重 /t 数量 / 辆 
小 型 车 2 40 
中 型 车 4 10 
大 型 车 8 29 
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矶 排放 因 子 "， 得 到 小 型 车 、 中 型 车 和 大 型 车 
的 CO, 排放 因子 分 布 为 326.88、501.64 和 
514.03 g/km. 

因此 ， 不 同 车 型 的 营运 车 辆 的 碳 排放 成 本 C, 
计算 如 公式 (4) 所 示 。 

© = DD E” L; (4) 


m=1i1=0j= 


KP, xt NEE; 若 车 型 为 m 和 车 辆 k 从 
ime, Wat = 1; 否则 xy = 0。E” 为 车 型 m 
的 碳 排放 因子 ，g/km。 工 ;为 任意 客户 点 i 到 j 的 距 
离 ，km。 

2.2.3” 货 损 成 本 计算 

本 研究 货 损 成 本 分 为 在 途 货 损 成 本 和 印 货 货 
损 成 本 

(1) 在 途 货 扣 成 本 。 指 运输 车 辆 向 客户 点 行 
驶 过 程 中 ， 货 物 因为 温度 、 湿 度 等 原因 造成 损失 
而 产生 的 成 本 。 生 鲜 农产品 的 特征 决定 了 它们 的 
货 损 成 本 和 时 间 成 正比 ， 所 以 假设 运输 过 程 中 单 
位 时 间 和 重量 货物 的 货 损 系数 : 

6, =(1- et) (5) 

其 中 ，s 是 产品 的 腐败 率 ，%， 值 越 大 ， 单 
人 T pn 
最 早 开始 服务 时 间 点 ，h; TY ,为 客户 点 i 前 一 个 
点 的 最 晚 服务 时 间 点 ，h。 

在 途 货 损 成 本 Ci 的 计算 公式 为 : 


= VY Sx p00, (6) 

AP, x HEE: 车 客户 i 点 由 mm 车 型 
的 车 辆 服务 ， 则 x” = 1; 和 否则 xz = 0。p 为 生 鲜 
货物 的 单价 ， 元 。0; 为 配送 车 辆 抵达 i 点 时 的 载 
WE, %o 

(2) BUGLE HAMAS. AEREA A 
后 ， 因 为 打开 车 厢 门 搬运 货物 等 原因 造成 货物 损 
耗 而 产生 的 成 本 。 假 设 御 货 时 单位 时 间 重 量 的 货 
物 的 货 损 系 数 为 : 

6,=(1- e™) (7) 

其 中 ，7, 为 客户 点 i 的 服务 时 间 ，h。 

印 货 时 的 货 损 成 本 CP 为 公式 (8)， 总 货 损 
成 本 C; 计 算 为 公式 (9): 


C? = > Sxp0.(0, 一 q) (8) 

C, = Ci + C? (9) 

其 中 ，g, 为 客户 点 i 的 货物 需求 量 ，kg。 
2.24 车 辆 固定 成 本 计算 

车 辆 的 固定 成 本 主要 包括 汽车 的 使 用 费用 、 
行驶 和 外 货 时 的 人 力 成 本 以 及 配送 车 辆 的 固定 局 
用 成 本 。 和 车 辆 的 固定 成 本 C, 公 式 为 : 


M K n n 
CS ES Dn ty (us + Pa) + 
M K n 
LEDs (a+r) $SS e 
m=1k=0i=0 


(10) 

其 中 ，M 为 车 辆 大 的 单位 时 间 使 用 费用 ， 元 ; 

P, 为 使 用 m 类 型 车 辆 的 单位 时 间 人 力 成 本 ， 元; 
6 为 车 型 m 的 固定 发 车 费用 ， 元 ; XW 为 0~1 变 
E: 若 车 型 为 m 的 车 辆 从 配送 中 心 到 j 点 ， 则 
XW 1， 否则 X% = 0; YY 为 0~1 变 量 : FEM 
为 m 的 车 辆 经 过 i 点 ， 则 Y =1, B Y? = 0。 


2.3 配送 路 径 多 目标 优化 模型 建立 


在 考虑 行驶 速度 、 车 辆 碳 排放 和 货 损 成 本 、 
时 间 窗 、 多 车 型 等 情况 下 ， 满 足 所 有 顾客 的 需 
求 ， 以 行驶 距离 C,、 碳 排放 成 本 C,、 货 损 成 本 
C; 和 车 辆 固定 成 本 C, 最 低 为 目标 ,综合 考虑 上 
述 四 个 目标 函数 ， 生 鲜 电 商 配 送 路 径 优 化 的 多 目 
标 联合 优化 模型 如 下 所 示 : 


minC, sS Ko (11) 
i=0j=0k 

minC, = Sye “BP + L; (12) 
m=1i=0j7=0 


min C, = C; + C$ = $ Sepo, (TTE) 
m=1i=0 
M n 
Dx" po:(O, - 4) Pu (13) 
=0 
mney -2 


m=1k= 

M K n 

DDD res (GHP) * yy eee, 
m=1k=0i=0 

(14) 


m=1k= 
其 中 ，xw 为 0~1 变 量 : EKA BM] 


ti (u, z P„) + 


ChinaXiv 合 作 期 刊 


143 


202302.00194v1 


chinaXiv 


点 ， MEJ =1, 否则 xw =0. 
3 改进 蚁 群 算法 设计 
3.1 改进 蚁 群 算法 设计 的 基本 思想 


蚁 群 算法 最 初 用 于 求 旅行 商 问题 (Traveling 
Salesman Problem, TSP) 等 单 目标 优化 问题 “， 
其 良好 的 鲁 棒 性 的 、 正 反馈 性 °°" 和 分 布 式 特 点 
可 较 好 地 解决 配送 路 径 优化 问题 S 。 本 研究 通 
过 采用 主要 目标 法 和 新 的 编码 方式 ， 将 该 问题 转 
换 为 可 求解 的 单 目标 问题 ; 接着 在 蚁 群 算法 中 融 
合 协同 进化 机 制 、 自 适应 信息 素 更 新 策略 和 局 部 
搜索 等 机 制 来 提高 算法 的 求解 效果 ; 最 后 设计 合 
适 的 适应 度 计算 方法 和 停 小 规避 策略 ， 增 强 该 算 
法 跳出 局 部 最 优 的 能 
3.1.1 主要 目标 法 

多 目标 模型 的 求解 通常 有 两 种 : 第 一 种 是 通 
过 线性 加 权 法 或 主要 目标 法 将 多 目标 优化 问题 转 
换 为 单 目 标 问 题 ; 第 二 种 是 通过 Pareto 方法 进行 
多 目标 同时 优化 ， 得 到 一 个 最 优 解 集 。 行 驶 距离 
最 短 、 生 鲜 农 产品 新 鲜 度 和 碳 排放 量 存 在 那 反 关 
系 7 。 对 于 冷 链 企业 来 说 ， 针 对 不 同 的 订单 情况 
需要 采取 不 同 的 优化 目标 ， 以 实现 成 本 最 小 和 有 僵 
利 最 大 。 因 此 ， 本 研究 采用 主要 目标 法 ， 分 别 以 
配送 总 成 本 最 小 、 货 损 成 本 最 小 、 车 辆 碳 排 放 和 
车 辆 固定 成 本 最 小 之 一 作为 主要 目标 ， 其 余 目 标 
满足 企业 的 要 求 即 可 。 基 于 上 述 分 析 ， 将 多 目标 
优化 问题 转化 为 单 目标 优化 问题 ， 当 行驶 距离 为 
主要 优化 目标 时 : 

min Cist. {C<scC<x0o,i=2,4} (15) 

其 中 ，56,、 忆 是 生 鲜 农产品 的 货 损 成 本 的 最 
高 和 最 低 约 束 ; 6,、 忆 是 车 辆 配送 过 程 中 的 碳 排 
放量 的 最 高 和 最 低 约束 ; 6 已 是 车 辆 模型 对 配 
送 过 程 中 ， 车 辆 固定 成 本 的 最 高 和 最 低 约 束 。 
3.1.2 ”种群 编码 和 解码 

在 利用 蚁 群 算法 求解 物流 路 径 优 化 问题 时 ， 


自然 数 编码 方式 。 改 进 凡 和 群 算法 编码 采用 如 下 设 
th: 假设 一 个 配送 中 心 为 2 个 客户 点 提供 配送 ， 
其 中 配送 中 心 有 M 种 类 型 车 辆 供 选择 ， 每 类 车 辆 
数目 为 K 辆 。 采 用 nn 维 实数 向 量 表示 种 群 中 每 个 
个 体 的 状态 ,使 用 mm 维 实数 向 量 来 表示 配送 车 辆 
类 型 ， 因 此 ， 当 进行 解码 操作 时 ， 首 先 按 顺序 从 
m 维 问 量 中 取得 车 辆 型 号 ， 接 着 根据 该 型 号 车 辆 
的 载重 确定 该 车 辆 需要 服务 的 客户 点 。 

车 辆 选择 和 客户 分 配 步 又 如 下 : 

(1) 得 到 初始 路 径 ; 

(2) 选取 最 大 承载 重量 车 辆 执行 配送 任务 ; 

(3) 判断 当前 车 辆 剩余 容量 能 否 满足 当前 客 
户 点 需求 ， 若 能 ， 转 步骤 (4), AMR 
步骤 (5); 

(4) 将 当前 客户 点 分 配给 当前 车 辆 ， 若 仍 有 
客户 未 服务 ， 转 步骤 (3)， 和 否则 输出 路 径 分 配 结 
果 ， 算 法 结束 ; 

(5) 选取 次 容量 作为 配送 车 辆 ， 转 
步骤 (3). 

3.1.3 ”种 群 协同 进化 机 制 

为 进一步 提高 算法 的 求解 效果 ， 本 研究 引入 
WERA (Ant Colony System, ACS) 和 最 大 最 
/\\ WEE AA (Max-Min Ant System, MMAS ) 种 
群 交 流 机 制 所 。 协 同 进 化 机 制 可 以 使 两 个 种 群 
从 算法 开始 到 结束 完全 联系 在 一 起 ， 不 再 单独 地 
进行 解 的 生成 和 迭代 ， 而 是 形成 两 类 蚂 蚊 种 群 相 
互 学 习 和 共同 促进 的 局 面 ， 使 得 算法 能 以 较 少 的 
迭代 次 数 找到 最 优 解 。 

对 于 协同 种 群 的 信息 素 交 流 ， 当 采取 固定 的 
交流 间隔 时 ， 会 导致 种 群 多 样 性 下 降 ， 产 生 的 结 
果 和 一 个 种 群 类 似 ; 当 信息 交流 间隔 过 长 时 ， 会 
减少 两 类 蚁 群 协 同 进化 的 效果 。 因 此 ， 本 研究 采 
用 基于 相似 度 的 交流 频率 策略 中， 对 交流 频率 
进行 自 适应 调节 ， 具 体 如 下 。 

当 两 类 蚁 群 算法 完成 本 次 迭代 ， 并 更 新 自身 


确定 解 的 表示 方法 是 一 项 非常 关键 的 工作 。 目 前 
普遍 使 用 的 编码 方式 有 二 进 制 编码 和 自然 数 编 
码 。 由 于 本 人 研究 涉及 到 多 车 型 的 表示 问题 ,采用 


计 息 素 后 ， 首 先 将 两 类 蚁 群 算法 的 全 局 最 优 路 径 
取 交 集 ， 得 到 两 类 蚁 群 算法 的 相似 度 。 当 相似 度 
较 低 时 ， 两 类 蚁 群 的 多 样 性 较 好 ， 不 必 进 行 交 
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流 。 当 相似 度 过 高 时 ， 两 类 蚁 群 路 径 趋 于 一 致 ， 
多 样 性 较 差 ,陷入 局 部 最 优 ， 应 进行 信息 素 交 
流 ， 将 两 类 蚁 群 的 信息 素 进行 交换 ”。 根 据 公 
式 (16) 判断 是 否 进 行 交流 。 


1S8,<S 
L= AB (16) 
| else 


其 中 ， 工 为 判断 是 否 进行 信息 素 交 流 的 逻辑 
fils S, 是 两 类 种 群 最 优 解 的 相似 度 ; 5 是 进行 交 
PLR. 

3.14 公共 路 径 奖 励 机 制 

公共 路 径 可 以 看 作 是 全 局 最 优 解 的 组 成 部 
分 ， 因 此 可 以 充分 利用 这 个 因素 ， 在 算法 前 期 ， 
增加 奖励 的 信息 素 ， 达 到 加 快 收敛 速度 的 目的 ; 
在 算法 后 期 ， 随 着 迭代 次 数 的 增多 ， 增 加 的 信息 
素 会 越 来 越 少 ， 从 而 增加 算法 后 期 的 多 样 性 。 

当 两 类 种 群 进行 交流 时 ， 对 种 群 的 公共 路 径 
的 信息 素 按照 公式 (17) 进行 更 新 : 

(t+ DPD) + Ay + (FT) (17) 


KP, 为 客户 点 的 总 数量 ， 个 ; iter 为 当 
前 迭代 次 数 ， 次 ; p 为 蚁 群 信息 素 的 蒸发 率 。 
3.1.5 ”局 部 搜索 (2-opt) 机 制 

本 研究 中 采用 该 算 子 来 加 快 收敛 速度 。 因 为 
在 不 同 车 辆 的 配送 路 径 间 进行 2-opt 优化 时 ， 
可 能 会 发 生 超载 现象 。 因 此 本 研究 采用 路 径 内 
2-opt 优 化 方式 。 优 化 方法 如 图 1 所 示 。 完 成 
每 次 操作 后 ， 计 算 该 条 配送 路 径 是 否 变 短 ， 若 路 
径 得 到 改进 ， 则 更 新 路 径 。 如 此 重复 ， 直 到 对 该 
条 线路 内 所 有 客户 都 进行 2-opt 操 作 。 


“i-l j} i) ”» i-l i | 


jt j” jtly <j) 
图 1 2-opt 优 化 示意 图 


Fig. 1 Schematic diagram of 2-opt optimization 


3.1.6 ”适应 度 计算 

适应 度 是 个 体 在 种 群生 存 的 优势 程度 度量 ， 
使 用 适应 度 函 数 来 进行 计算 ， 在 个 体 进 化 过 程 中 
起 着 重要 的 推动 作用 ， 适 应 度 太 大 会 导致 算法 过 
早 收 僵 ， 太 小 则 会 使 求解 过 程 缓慢 '， 因 此 ， 


选择 合适 的 适应 度 值 很 重要 。 在 求解 车 辆 路 径 优 
化 问题 时 ， 目 的 是 使 得 总 目标 函数 最 小 且 非 负 ， 
因此 为 了 使 所 求 优化 目标 函数 最 小 ， 则 应 使 适应 
度 值 最 大 。 

在 蚁 群 算法 的 求解 中 ,个 体 被 保留 的 可 能 性 
与 适应 度 值 呈正 相关 ， 适 应 度 值 较 大 的 个 体 越 有 
可 能 被 保存 下 来 ， 所 以 适应 度 忆 和 模型 总 成 本 的 
关系 应 为 倒数 。 本 人 研究 以 C, 代 表 第 i 条 路 径 所 产 
生 的 总 成 本 ,适应 度 值 计算 如 公式 (18) '”。 

1 

er (18) 
3.1.7 ”停滞 规避 策略 

基本 蚁 群 算法 在 迭代 一 定 次 数 后 ， 可 能 会 陷 
人 局 部 最 优 ， 并 在 较 差 路 径 上 积累 较 多 信息 素 ， 
进而 出 现 停滞 现象 ”%。 但 基本 蚁 群 算法 本 身 缺 
乏 跳 出 局 部 最 优 的 能 力 ， 因 此 本 研究 针对 该 问 
题 ， 采 取 了 一 种 停滞 检测 和 跳出 策略 ? 。 

算法 停 沾 检测 办 法 是 计算 更 优 解 的 次 数 。 妆 
在 一 定 和 迭代 次 数 内 ， 仍 未 出 现 较 优 解 ， 则 判定 算 
法 陷 人 局 部 最 优 ， 有 具体 检测 判断 公式 如 下 : 

_ n count > T 


(19) 
0 else 


其 中 , 工 是 判断 算法 是 否 停滞 的 0~1 布 尔 
值 ; count 是 判断 算法 停 沾 的 浆 值 。 

针对 算法 停滞 的 具体 措施 如 下 : 当 检 测 到 算 
法 中 全 局 最 优 值 陷入 停滞 时 ， 首 先 将 两 类 蚊 群 中 
所 有 路 径 上 的 信息 素 置 为 初始 值 ， 接 着 科 选 出 两 
类 疏 群 停 沛 前 所 产生 最 优 路 径 的 公共 路 径 ， 最 后 
奖励 公共 路 径 上 的 信息 素 ， 公 式 如 下 : 


T(t +1)= T +AT; (20) 
1 
— if (i,j BestTour 

wenu (21) 
0 else 


其 中 ，zm 为 每 条 路 径 上 的 起 始 信息 素 浓度 ; 
At; 是 两 类 蚁 群 分 别 由 公式 (21) 得 出 的 信息 素 
增 量 。 强 化 公共 路 径 的 信息 素 会 增强 算法 停滞 后 
再 次 运行 的 导向 性 ， 并 可 以 使 算法 的 搜索 侧重 点 
放 在 非 公 共 路 径 的 城市 序列 中 ， 使 得 算法 易于 发 
现 更 优 的 路 径 ， 保 证 算法 在 后 续 的 循环 中 能 有 效 
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跳出 局 部 最 优 。 
3.2 改进 蚁 群 算法 步骤 

改进 蚁 群 算法 的 流程 见 图 2。 图 2 具体 步 又 
如 下 : 

(1) 参数 初始 化 。 

(2) 将 种 群 划分 为 两 组 ， 分 别 按照 最 大 最 小 
疏 群 算法 和 普通 蚁 群 算法 更 新 信息 素 ， 求 解 种 群 
的 适应 度 。 

(3) 进入 协同 交流 阶段 ， 计 算 两 类 蚁 群 的 路 
径 相 似 度 ， 接 着 根据 公式 (17) 判断 是 否 进行 交 
流 ， 若 是 ， 转 步骤 (4)， 否 则 转 步 又 (5)。 

(4) 进入 种 群 交流 阶段 ， 找 出 种 群 的 公共 路 
径 ， 按 照 公式 (20) 奖励 公共 路 径 的 信息 素 。 

(5) 根据 公式 (19) 对 算法 是 否 陷 入 局 部 最 
优 进 行 检查 ， 若 是 ， 转 步骤 (6) ,否则 转 
步骤 (7)。 

(6) 根据 公式 (20) 对 算法 陷入 局 部 最 优 做 
出 相应 处 理 。 

(7) 交换 信息 素 ， 达 到 跳出 局 部 最 优 的 
效果 。 

(8) 判断 当前 迭代 次 数 是 否 大 于 或 等 于 最 大 
迭代 次 数 ， 若 是 ， 输 出 最 优 解 ， 算 法 结束 ， 否 则 
FEBRILE GE (2)。 


4 数值 模拟 


为 验证 本 研究 提出 的 生 鲜 农产品 配送 路 径 优 
化 模型 和 改进 的 蚁 群 算法 ， 将 带 时 间 窗 的 车 辆 路 
径 问题 (Vehicle Routing Problems With Time 
Windows, VRPTW) 的 Solomon 数据 集 算 例 六 
进行 改造 ， 使 用 Matlab2010 进 行 仿真 实验 。 设 有 
一 个 生 鲜 农产品 配送 中 心 和 100 个 客户 ， 各 客户 
点 服务 时 间 、 时 间 窗 如 原本 算 例 3 所 述 ， 配 送 
中 心 的 冷 链 配 送 车 辆 的 载重 量 分 别 为 2、4 和 8t， 
各 类 型 车 辆 数量 分 别 为 5、3 和 2 辆 。 假 设 配 送 货 
车 以 35 km/h 的 速度 匀速 行驶 。 


4.1 参数 设置 


改进 模型 和 算法 实验 以 Solomon 数据 库 的 


始 


y 
参数 初始 化 
路 径 构建 
ozs 


¥ x 
根据 ACS 更 新 根据 MMAS 更 新 
Y Y 
信息 素 更 新 信息 素 更 新 
| | X 
| 协同 交流 
EE E 
停滞 规避 
b d 
NC<=NCMAX 


l 党 


| 输出 最 优 路 径 | 
CE 
图 2 改进 蚁 群 算法 流程 图 
Fig. 2 Flow chart of improved ant colony algorithm 
C105, R105 和 RC105 为 例 ， 算 例 C105 属于 集中 
分 布 ，R105 的 客户 坐标 属于 随机 分 布 ，RC105 
属于 混合 分 布 。 参 考 其 他 人 研究 结果 “””， 并 经 过 
大 量 实验 后 ， 参 数 设置 如 下 ; 令 客 户 需求 量 的 
个 重量 单位 为 10 kg， 最 终 总 配送 重量 为 57t， 更 
符合 实际 需求 ; 为 简化 计算 ， 此 次 数值 计算 假设 
生 鲜 农产品 类 型 单一 ， 其 中 产品 腐败 率 s=0.06， 
货物 单价 为 1 元 /kg。 计 算 固 定 成 本 时 ， 冷 链 配 
送 车 辆 的 人 力 成 本 、 启 用 成 本 和 管理 成 本 如 表 2 
所 示 ， 碳 税 计算 价格 为 C=0.05 元 /kg 。 通 过 反 
复 实 验 ， 改 进 蚁 群 算 法 的 参数 分 别 设 置 为 蚂 尽 数 
m 为 50 个 ,信息 素 重 要 程度 启发 因子 a=1， 启 发 
函数 重要 程度 因子 B=3， 信 息 素 挥发 因子 p 为 
0.85， 控 制 转移 规则 参数 m 为 0.5。 迭 代 次 数 设置 
为 100 次 。 
表 2 车 辆 的 人 力 、 启 用 和 管理 成 本 


Table 2 Labor, start-up and management costs 


of vehicles 
车 辆 类 型 小 型 车 中 型 车 ” 大 型 车 
AFI RAS (FE + km) 20 30 40 
启用 成 本 /( 元 :次 ”) 300 400 500 
管理 成 本 /( 元 :次 ”) 60 80 100 
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4.2 算法 性 能 对 比 
4.2.1 不 同 优化 目标 结果 比较 

为 验证 模型 的 有 效 性 ， 以 本 研究 提出 的 模型 
为 基础 ， 采 用 改进 后 的 Solomon 数据 集中 的 
C105 算 例 ， 使 用 改进 蚁 群 算法 求解 模型 ， 分 别 
以 总 路 径 最 短 ， 产 品 新 鲜 度 最 大 ， 碳 排放 以 及 车 
辆 固定 成 本 最 小 为 主要 目标 ， 其 他 目标 转 为 约束 
进行 求解 ， 具 体 结果 如 表 3 所 示 。 

表 3 不 同 优化 目标 的 结果 比较 


Table 3 Comparison of results of different optimization 


objectives 
Re Nay wu 车 辆 固 
主要 优化 目标 MEARKE WHEE 产品 新 鲜 度 、 、 
定 成 本 
总 路 径 /km 937.94 968.91 1105.79 965.02 
碳 排 放 /元 5136.02 4961.48 6150.46 5095.24 


产品 货 损 成 本 /元 6126.35 5986.93 4081.78 5875.46 


车 辆 固定 成 本 /元 7739.21 7573.20 8219.79 7500.87 


从 表 3 中 可 以 看 出 ， 当 选择 不 同 的 主要 优化 
目标 时 ， 该 模型 均 能 在 该 目标 求 得 最 优 解 : 尤其 
对 冷 链 配送 的 关键 指标 产品 新 鲜 度 有 了 很 大 提 
升 ， 最 高 降幅 达到 36%; 以 路 径 总 长 度 为 优化 目 
标 时 ， 取 得 的 结果 接近 算 例 最 优 解 ， 以 多 车 型 的 


碳 排 放 计 算 为 目标 时 ， 最 多 减少 了 19% 的 碳 排 
放 ; 对 于 车 辆 固定 成 本 ， 降 幅 最 高 达到 了 5%。 
同时 几 组 解 的 结果 证 明 客户 规模 较 大 情况 下 ， 车 
辆 固定 成 本 、 生 鲜 农产品 新 鲜 度 和 碳 排放 量 存在 
悖 反 关 系 ， 多 个 优化 目标 很 难 同时 达到 最 优 。 例 
如 ， 寿 满 足 其 产品 货 损 成 本 最 低 ， 与 最 优 值 相 
比 ， 使 得 碳 排放 增幅 23%， 车 辆 固定 成 本 增加 
10%; 从 配送 成 本 的 构成 比例 来 看 ， 占 比 最 高 的 
车 辆 固定 成 本 平均 比例 达到 39%， 其 次 是 产品 新 
鲜 度 平均 占 比 27%。 因 此 ， 上 述 数 值 计算 表明 ， 
改进 模型 可 以 在 对 应 的 优化 目标 上 取得 最 优 结 
R, 该 模型 有 效 ; 对 于 冷 链 电 商 企业 ， 在 面 对 大 
规模 订单 时 很 难 同时 在 多 个 优化 目标 取得 最 优 ， 
决策 者 应 根据 自身 战略 的 需要 和 市 场 与 客户 不 同 
要 求 ， 选 择 特定 主要 目标 进行 优化 ， 进 而 得 出 适 
合 的 配送 方案 。 
4.2.2 不 同 算法 的 结果 比较 

将 本 研究 提出 的 算法 与 基本 蚁 群 算法 进行 对 
比 实验 ， 算 例 选 取 改 编 后 的 C105 算 例 ， 实 验 参 
数 不 变 ， 且 分 别 以 改进 蚁 群 算法 和 基本 蚁 群 算法 
对 其 进行 求解 ， 分 别 以 路 径 总 长 度 、 碳 排放 、 产 
品 新 鲜 度 和 车 辆 固定 成 本 为 主要 优化 目标 ， 结 
如 表 4 所 示 。 


表 4 不 同 算法 求解 算 例 计算 结果 


Table 4 Calculation results of examples solved by different algorithms 


Ree B 路 径 总 长 度 碳 排放 (调整 ) 产品 新 鲜 度 车 辆 固定 成 本 
ue 改进 ACS 基本 ACS 改进 ACS 基本 ACS 改进 ACS 基本 ACS 改进 ACS 基本 ACS 
路 径 总 长 度 /km 937.94 1025.02 968.91 1051.46 1105.79 1127.91 965.02 1066.76 
碳 排 放 / 元 5136.02 6115.47 4961.48 6051.10 6150.46 6670.16 5095.24 6307.41 
产品 货 损 成 本 /元 6126.35 7573.59 5986.93 7984.48 4081.78 6597.01 5875.46 6824.19 
车 辆 固定 成 本 /元 7739.21 8380.63 7573.20 8353.51 8219.79 8433.54 7500.87 8079.81 


以 路 径 长 度 为 主要 优化 目标 时 ， 碳 排放 和 货 
损 成 本 降幅 最 大 ， 达 到 16% 和 19% ， 总 配送 成 本 
降低 13%; 以 碳 排放 为 主要 优化 目标 时 ， 碳 排放 
和 货 损 成 本 在 总 配送 成 本 减少 量 中 贡献 最 多 ， 所 
占 比 例 分 别 为 50% 和 27%， 总 配送 成 本 降幅 
16%; 以 货 损 成 本 和 车 辆 固定 成 本 为 主要 优化 目 
标 时 ， 总 配送 成 本 分 别 降低 14% 和 12%; 不 同 优 


化 目标 的 总 配送 成 本 平均 降幅 达到 14%， 改 进 蚁 
群 算法 结果 更 优 。 综 上 所 述 ， 本 人 研究 的 改进 蚁 群 
算法 可 以 降低 冷 链 配 送 成 本 。 
4.2.3 不 同 客户 分 布 算 例 结果 比较 

为 进一步 验证 本 研究 生 鲜 农产品 配送 路 径 优 
化 模型 的 有 效 性 ， 采 用 多 类 型 算 例 C105、R105 
和 RC105 分 别 进行 试验 ， 三 种 客户 分 布 分 配 的 车 
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辆 数量 和 类 型 一 致 ， 其 他 试验 参数 也 一 致 ， 以 路 
径 总 长 度 为 主要 优化 目标 ， 试 验 结果 如 表 5 所 
示 ， 三 种 分 布 路 线 如 图 3~ 图 5 所 示 。 

表 5 不 同 客户 分 布 算 例 计算 结果 


Table 5 Calculation results of different customer distribu- 


tion examples 


分 布 类 型 路 径 总 WAE ”产品 货 损 车 辆 固定 ”总 成 
长 度 /km ğum ”成 本 /元 ”成 本 /元 ”本 /元 


C105 937.94 5136.02 6126.35 7739.21 19,939.53 
R105 1129.97 6659.98 7633.96 8671.55 24,095.47 
RC105 1185.22 6928.16 7615.82 8668.36 24,397.58 


距离 /km 
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图 3 C105 算 例 配送 路 径 示意 图 


Fig. 3 C105 example distribution path diagram 


距离 /km 


图 4 R105 算 例 配送 路 径 示意 


Fig. 4 R105 example distribution path diagram 

由 表 5 可 以 看 出 ，C 类 算 例 的 配送 距离 、 和 车 
辆 固定 成 本 、 货 损 成 本 和 碳 排 放量 在 所 有 类 型 中 
最 低 。 主 要 由 于 C 类 客户 分 布点 主要 集中 分 布 在 
儿 个 区 域 , 行驶 距离 较 短 ， 因 此 碳 排 放量 相对 较 
低 ， 并 且 货 物 中 途 货 损 成 本 保持 在 一 个 较 低 值 。 
但 R 类 和 RC 类 总 配送 距离 、 车 辆 固定 成 本 、 货 


AE X 
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图 5 RC105 算 例 配 送 路 径 示意 图 
Fig. 5 RC105 example distribution path diagram 
损 成 本 等 均 大 于 C 类 ， 主 要 由 于 这 两 种 类 型 客户 
均 随 机 分 布 ， 使 得 车 辆 配送 距离 增 大 ， 而 且 总 行 
驶 时 间 较 长 ， 总 配送 成 本 也 随 着 分 布 复杂 度 的 增 
加 而 增加 。 从 最 终结 果 对 比 来 看 ， 所 有 类 型 的 成 
本 里 面 ， 碳 排放 和 货 损 成 本 占 比 相对 较 高 ， 两 类 
成 本 占 比 均 达到 68% 以 上 。 碳 排放 随 着 分 布 的 复 
杂 化 涨幅 最 高 ， 平 均 涨 幅 16.5%。 说 明 无 论 客户 
的 分 散 程度 如 何 ， 车 辆 固定 成 本 和 产品 货 损 成 本 
应 是 企业 主要 关注 的 目标 。 因 此 要 降低 总 物流 成 
本 ,最 主要 就 是 合理 规划 配送 策略 ， 适 当 增 加 配 
送 中 心 ,减少 客户 的 分 散 度 。 合 理 降低 生 鲜 农 产 


品 货 损 。 
4.3 灵敏 度 分 析 


为 了 分 析 腐 败 率 和 车 辆 速度 对 冷 链 配送 成 本 
的 影响 ， 选 择 货 损 成 本 作为 主要 优化 目标 ， 在 保 
持 其 他 条 件 不 变 的 前 担 下 ， 令 各 参数 以 倍数 递增 
( 减 ) 两 次 ， 记 录 每 次 的 主要 优化 目标 值 并 分 析 。 
4.3.1 腐败 率 对 成 本 的 影响 

腐败 率 变 化 对 结果 的 影响 如 图 6 所 示 。 当 腐 
败 率 变化 时 ， 对 产品 的 货 损 成 本 影响 最 明显 ; 腐 
败 率 越 高 ， 产 品 的 货 损 成 本 越 高 ， 当 腐败 率 每 增 
长 10% 时 ， 货 损 成 本 增加 11%。 虽 然 腐 败 率 和 配 
送 总 成 本 无 明显 的 比例 关系 ,但 腐败 率 每 增 
( 减 ) 10% 时， 配送 总 成 本 均 增 加 。 因 此 ， 货 损 
成 本 并 不 是 越 低 越 好 ， 对 于 企业 来 说 ， 应 以 货 损 
成 本 为 主要 优化 目标 ， 使 货 损 成 本 保持 在 一 个 合 
理 的 水 平 ， 配 送 总 成 本 最 低 。 
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图 6 腐败 率 变化 对 货 损 和 总 成 本 的 影响 
Fig. 6 Impact of speed change on cargo damage cost and 
total cost 

4.3.2 ”车 辆 速度 变化 对 成 本 的 影响 

车 辆 速度 增 减 对 优化 目标 的 影响 如 图 7 所 
示 。 随 着 速度 的 增加 ， 配 送 车 辆 的 在 途 时 间 缩 
短 ， 使 得 碳 排放 减少 ， 最 终 总 成 本 呈 下 降 趋势 ; 
碳 排 放 与 速度 增 减 有 明显 的 比例 关系 ， 随 着 速度 
增加 而 降低 。 但 在 速度 增加 后 期 总 成 本 出 现 边际 
效应 ， 总 配送 成 本 减少 放 缓 。 因 此 ， 冷 链 企业 可 
以 通过 适当 提高 速度 来 减少 碳 排放 量 和 配送 
成 本 。 
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一 去 一 总 配送 成 本 ”一 里 一 碳 排 放 成 本 
图 7 速度 变化 对 碳 排放 和 总 成 本 的 影响 
Fig. 7 Impact of speed change on carbon emission cost and 


total cost 
5 结论 与 展望 


本 研究 以 VRPTW 模型 为 基础 ， 提 出 了 基于 
载重 的 多 车 型 碳 排 放 和 车 辆 速度 计算 模型 ， 建 立 
带 有 车 辆 启动 固定 成 本 、 路 径 长 度 和 货 损 成 本 的 
生 鲜 农产品 配送 模型 ， 并 设计 了 基于 协同 进化 思 
想 的 改进 蚁 群 算法 。 数 值 模 拟 结果 证 明 : 


(1) 相 比 单车 型 配送 模式 ， 本 研究 改进 后 的 
模型 可 以 有 效 降低 碳 排放 成 本 和 配送 总 成 本 。 根 
据 生 鲜 农产品 的 腐败 率 、 客 户 重 要 程度 、 配 送 场 
景 的 不 同 ， 采 用 相应 的 优化 目标 ， 可 以 使 企业 、 
客户 的 利益 最 大 化 。 

(2) 改进 的 蚁 群 算法 与 基本 蚁 群 算 法 相 比 ， 
可 以 有 效 求解 大 规模 客户 算 例 ， 配 送 总 成 本 平均 
降幅 达 14%， 进 一 步 降 低 企业 配送 成 本 。 

(3) 对 于 不 同 分 布 类 型 的 客户 订单 ， 客 户 集 
中 分 布 时 的 配送 成 本 最 低 ， 部 分 集中 时 需要 的 成 
本 最 高 ， 企 业 应 通过 增加 配送 中 心 等 方法 来 降低 
成 本 ， 提 高 配送 效率 。 

(4) 对 于 腐败 率 不 同 的 生 鲜 农产品 ， 货 损 成 
本 最 低 并 不 能 保证 总 配送 成 本 最 低 ， 应 使 货 损 成 
本 保持 在 一 个 合适 的 区 间 。 

(5) 车 辆 速度 越 高 ， 配 送 总 成 本 不 一 定 越 
低 ， 应 选择 合适 的 速度 降低 总 成 本 。 

本 研究 通过 提供 了 多 种 优化 目标 的 选择 ， 使 
得 企业 可 以 根据 自身 需要 和 发 展 战 略 ， 权 衡 不 同 
优化 目标 的 利 次 ， 得 出 适合 的 配送 方案 ; 同时 为 
智慧 农业 供应 链 的 配送 路 径 规划 体系 提供 了 新 的 
理论 和 数据 支撑 ， 有 较 强 的 应 用 价值 和 理论 
价值 。 

本 研究 仅 考虑 了 天 气 较 好 情况 下 的 车 辆 速度 
变化 ， 没 有 考虑 其 他 天 气 情况 下 对 配送 效率 的 影 
响 。 同 时 只 研究 了 单 温度 和 单 配 送 中 心 下 的 冷 链 
配送 车 辆 路 径 问 题 ， 没 有 研究 多 温 共 配 和 多 配送 
中 心 的 农产品 配送 问题 ”5;， 且 仅仅 使 用 主要 目 
标 法 求解 多 目标 模型 ， 没 有 考虑 其 他 方法 。 因 
此 ， 以 上 问题 将 是 下 一 步 的 研究 内 容 。 
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Time-Varying Heterotypic-Vehicle Cold Chain Logistics 
Distribution Path Optimization Model 


LIU Siyuan'*, CHEN Tian'en”, CHEN Dong’, ZHANG Chi’, WANG Cong’ 


(1. School of Computer, Electronics & Information, Guangxi University, Nanning 530004, China; 
2. National Engineering Research Center for Information Technology in Agriculture(NERCITA), Beijing 10097, China) 


Abstract: In view of the problems of constant speed and single carbon emission calculation method in the distribution model of 
fresh agricultural products in the transportation link of agricultural supply chain, combined with the time-varying characteris- 
tics of road network and the new multi vehicle carbon emission calculation method, this study put forward the distribution route 
optimization model of fresh agricultural products with four optimization objectives, which were the distribution distance, multi 
vehicle carbon emission, goods loss and vehicle fixed cost. In this model, the calculation of fuel consumption and carbon emis- 
sion in the model would be affected by many factors, among which the load is the most important factor: Firstly, the average fu- 
el consumption per 100 km of different trucks was calculated, then the CO, emission factors of various trucks were calculated 
according to the carbon balance principle, and finally the average value of the results of each truck was taken as the carbon 
emission factor of the vehicle. According to those characteristics of the model, an improved double strategies co-evolutionary 
ant colony system (DC-ACS) was proposed. In this study, the main method was used to transform the problem into a solvable 
single objective problem. Then, the ant colony algorithm combined the coevolution mechanism, adaptive pheromone update 
strategy and local search mechanism were used to improve the solution effect of the algorithm. Finally, an appropriate fitness 
calculation method and stagnation avoidance strategy were designed to enhance the ability of the algorithm to jump out of local 
optimization. The C105 example of Solomon dataset was solved by using the improved ant colony algorithm. The optimal solu- 
tions on the four optimization objectives were 937.94 km, 4961.48 CNY, 4081.78 CNY and 7500.87 CNY respectively, which 
proved the effectiveness of the model proposed in this study. Based on the effectiveness of the model, the experimental results 
showed that the total distribution cost of the improved ant colony algorithm reduced by more than 14% on average compared 
with the basic ant colony algorithm on the four optimization objectives, which proved that the improved ant colony algorithm 
had more advantages. The improved ant colony algorithm was used to solve large-scale examples with different distributions: 
centralized, random and mixed. The optimal total costs were 19939.53 CNY, 24095 CNY and 24397.58 CNY, respectively. To 
sum up, the proposed model and algorithm could provide a good reference for the urban distribution path decision-making of 
cold chain logistics enterprises, a new idea to improve the distribution path optimization model and optimization method of 
smart agricultural supply chain, and a reference for enterprises to further expand their scale. 


Key words: cold chain logistics; path optimization; real-time information; ant colony optimization; Solomon dataset 
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